Resonant demagnetization of a dipolar BEC in a 3D optical lattice 
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We study dipolar relaxation of a chromium BEC loaded into a 3D optical lattice. We observe dipo- 
lar relaxation resonances when the magnetic energy released during the inelastic collision matches 
an excitation towards higher energy bands. A spectroscopy of these resonances for two orientations 
of the magnetic field provides a 3D band spectroscopy of the lattice. The narrowest resonance 
is registered for the lowest excitation energy. Its line-shape is sensitive to the on-site interaction 
energy. We use such sensitivity to probe number squeezing in a Mott insulator, and we reveal the 
production of three-body states with entangled spin and orbital degrees of freedom. 

PACS numbers: 03.75.Mn , 05.30.Jp, 67.85.-d 



Quantum dipolargascs have attracted much attention 
in recent years [lj-yl- Following seminal studies of Cr 
Bose Einstein Condensates (BEC s) [4| -|7|, and the recent 
production of Er Q and Dy 0, fief quantum gases, it 
is natural to study dipolar gases confined in optical lat- 
tices. Dipolar gases in lattices provide an ideal play- 
ground for studying quantum phase transitions in a sys- 
tem with long range interactions. The unique proper- 
ties of dipole-dipole interactions (DDIs) presents direct 
similarities with the Heisenberg model of quantum mag- 
netism [ll-14|, and lead to novel quantum phases dis- 
playing possible long range ordering in lattices 
The non-linear coupling between spin and orbital de- 
grees of freedom provided by dipolar interactions [iHI^ 
is particularly interesting. Although at large magnetic 
field, this coupling leads to fast dipolar relaxation losses 



for atoms in excited spin states |24(, we have recently 
demonstrated that working at low magnetic field enables 
the study of multicomponent quantum gases with free 
magnetization 0, (2(| . Here we extend this research to 
include dipolar particles trapped in optical lattices, and 
directly observe discrete, atom-number dependent cou- 
pling between spin and orbital degrees of freedom. 

We study the magnetization of 52 Cr atoms loaded in 
a deep 3D optical lattice. Under our experimental con- 
ditions, we produce a Mott insulator state with a core 
of two particles per site 27, 28j . We show that on-site 
dipolar relaxation is inhibited unless the released mag- 
netic energy matches a lattice band excitation, in which 
case we observe demagnetization resonances. The in- 
terplay between the anisotropics of dipolar interaction 
and of lattice sites leads to a resonance spectrum which 
depends on the magnetic field orientation. Measuring 
demagnetization as a function of the magnetic field for 
two orientations allows for spectroscopy of the 3D lattice 
band structure. The narrowest resonance is found to be 
sensitive to the on-site atom number distribution and re- 
veals the number squeezed distribution of the Mott state. 
Changing our experimental conditions, we prepare sites 
containing three atoms. Spin-orbit coupling DDI then 
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FIG. 1: (Color online) (a) Sketch: Dipolar relaxation only occurs 
provided the loss of Zecman energy AEzccman matches a gain in orbital 
energy (band excitation), (b) State composition at resonance after 6 
ms and after 30 ms of dipolar relaxation (absorption images after Stern- 
Gcrlach separation). 



produces three-body states that are expected to have en- 
tangled spin and orbital degrees of freedom. 

We produce a 52 Cr BEC (with spin S — 3) with typ- 
ically N = 1, 5 ■ 10 4 atoms by forced evaporation in the 
lowest energy state \ms = — 3) in a crossed optical dipolc 
trap derived from a fiber laser at 1075 nm [29[ (ms is the 
spin projection along the magnetic field). The BEC is 
loaded adiabatically in a 3D optical lattice, generated 
using 4 W of a single mode laser at A=532 nm. A 2D 
horizontal lattice is created by the interference of three 
beams: two beams are counter propagating, the third 
intersects the retro-reflected pair at 45°. The resulting 
2D lattice is rectangular with a periodicity of A/ (2 cos ^) 
along Ox direction and A/ (2 sin ^) along Oy direction. In 
the vertical direction an independent ID lattice is created 
by a vertical rctro-reflcctcd beam with horizontal polar- 
ization. The 1/e 2 beams radii at the BEC position are 
60(10) fim. The lattice depth is measured by Kapitza- 
Dirac diffraction ^ to be 25 E R (with E R = h 2 /2m\ 2 
the recoil energy), corresponding to first band excitation 
frequencies of uj x /2tt = 170(10) kHz, uj v /2tt = 55(5) kHz 
and lu z /2tt = 100(10) kHz. 

For most of the experiments presented here, we load 
the optical lattice adiabatically in 40 ms. Although the 
lattice geometry is complicated by its anisotropy, and we 
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FIG. 2: (Color online) Observation of resonances due to dipolar relaxation as a function of the magnetic field amplitude. The spectra (red 
bullets - after 30 ms of dipolar relaxation time, black diamonds - after 20 ms) are obtained respectively for a magnetic field along Ox and Oy, two 
orthogonal directions in the horizontal plane. The sketch on the right represents the two studied magnetic field orientations referred to a lattice 
site. The peaks represent the frequency positions and the calculated coupling amplitudes of the resonances leading to excitations of lattice bands 
labeled by numbers (n x , n y , n z ) . Shaded areas qualitatively represent energy broadening due to tunneling in excited bands. 



do not expect a fully adiabatic crossing towards the Mott 
state, this loading time scale allows to remain close to the 
ground state. This is confirmed by adiabatically ramp- 
ing down the lattice at the same speed, in which case a 
BEC with little excitations is recovered. For our experi- 
mental parameters, the predicted Mott state is made of 
a central region with double occupancy surrounded by a 
single occupancy shell. 

Once the lattice is loaded, atoms are adiabatically 
transferred from |mg = — 3} = |— 3} to the highest en- 
ergy state |3) by means of a 5 ms radio-frequency (rf) 
sweep (the efficiency of rf sweeps is 85 %). Once atoms 
are transferred to |3), dipolar relaxation is allowed. After 
a given time, the magnetization is again reversed using a 
second rf sweep (this step is performed only to optimize 
the detection), the optical lattice is switched off in I ms 
and the populations in the different magnetic states are 
monitored by a Stern Gerlach analysis (Fig. Q]b). 

We obtain two different behaviors as a function of the 
magnetic field amplitude B i.e. whether the Larmor fre- 
quency TiWLar = 9s^bB is smaller or larger than the lat- 
tice depth Vq (/zb is the Bohr magneton and g$ — 2 is the 
Lande factor). For large magnetic fields, TiWLar > Vo, we 
observe a non exponential loss dynamics: the atom num- 
ber rapidly decreases before reaching a non zero steady 
value and then no loss are observed for tens of ms. Losses 
are due to dipolar relaxation between atoms occupying 
the same site and they stop when only single occupancy 
sites remain. 

At low magnetic field, huitar <C Vo, no losses are ob- 
served, as atoms remain trapped after dipolar relaxation. 
We instead observe a demagnetization dynamics, that re- 
veals a transfer of atoms from |3) to |2) and then to lower 
energy magnetic states (Fig.[T]b). For a fixed relaxation 
time we registered the final population in |3) as a func- 
tion of the magnetic field amplitude. The spectra pre- 
sented in Fig. [2] show resonance peaks which we interpret 



as the consequence of total energy conservation: for a 
pair of atoms initially in the vibrational ground state of 
one lattice site, spin relaxation only occurs if the released 
magnetic energy matches a band excitation. We obtain 
two different spectra for two orthogonal directions of the 
magnetic field B. This shows that dipolar relaxation cou- 
ples the initial state to different lattice bands depending 
on the orientation of B relative to the geometry of a lat- 
tice site. 

To interpret our data, we consider a simple theoret- 
ical model with two particles in a purely 3D harmonic 
anisotropic potential (with frequencies uj x ,uj v ,uj z ), de- 
scribing one lattice site. Dynamics is determined by the 
dipole-dipole interaction Vdd(r) = ^ [r 2 Si • S 2 — 3(Si • 

r)(S2 • r)], where d 2 = /j,q (gs/is) 2 /47r, (j,q the magnetic 
constant; r = ri — r 2 is the relative coordinate; Sj is the 
spin operator for atom i. In presence of an harmonic con- 
finement, the center of mass and the relative motions are 
decoupled. DDIs affect only the relative motions, there- 
fore only the relative coordinates need to be considered, 
with eigenstate \n x ,n y ,n z ). For the magnetic state of the 
pair, we use the molecular basis \Sm>^m) where Sm is 
the total spin and tom is the spin projection. The energy 
is given by E n .„ lM = £\ ^(n, + \) + m M g s HBB where 
i = x, y, z. The particles are prepared in the initial state 
|i) = 0, 0, 0) ® |6,6). 

From this initial state, due to DDIs, three spin- 
relaxation channels a re p ossible with final spin states 
6,5), |6,4) and |4,4) [24|. The couplings to these three 
final states are described by the DDIs coupling operators 

Vl = 3SU 2 (x' + iy')z'/f 5 , V 2 = %^Sd 2 (x' + iy') 2 /f 5 

and V3 = ^^J~^Sd 2 (x' + iy') 2 /r 5 where the orientation 

of the magnetic field B is along the z'-axis 0| ( x',y' 
and z' are the relative coordinates). The first channel 
induces a loss of magnetic energy gs^BB, and the two 
others of 2gs/isB. 
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FIG. 3: (Color online) High resolution spectroscopy of the lowest 
frequency dipolar resonance: sensitivity to site atom number, (a) With 
two atoms per site (red bullets), the line shape is narrow, while with 
three atoms per sites (blue triangles), a multi-resonant structure reveals 
the production of three-body states with entangled spin and orbital de- 
grees of freedom (two peaks with a star). The solid lines (top) are fits, 
whereas the bottom dashed lines result from calculation (including an 
empirical broadening mimicking technical noise). Each data point is 
the average of five measurements, (b) Sketch to interpret the site occu- 
pancy sensitivity, arising from an on-site interaction shift depending on 
the orbital state in the lattice, (c) Resonance width (obtained by fit- 
ting the data with a single gaussian) as a function of the lattice loading 
time. 



As sketched in Fig. [T] the Zeeman energy is converted 
into excitation to a higher orbital state \n x , n y ,n z ). Due 
to the dependence of Vi,2.3 on x', y' and z', the selection 
rules and the coupling strengths to final orbital states 
strongly depend on the orientation of the magnetic field 
B. This explains the difference between the two regis- 
tered spectra (see Fig. [2]) for B along Ox and Oy respec- 
tively. DDIs couple the initial orbital state with orbital 
excited states with even total symmetry. Additional spe- 
cific selections rules are derived from the parity of V\ 2 3 
in the 3 directions x', y', z' . In practice, most of the reso- 
nances observed in Fig. [5] are attributed to channel 2 and 
3. The observed resonances corresponds to a rate of the 
order of 100 Hz in good agreement with a calculation 
of the coupling strengths (i\ Vi.,2,3 |/). The relatives cou- 
pling strengths are represented by the peaks amplitudes 
in Fig. [5] and we also give the their value in supplemental 
material. 

A striking observation in Fig. [2] is that resonances 
widen as the excitation energy increases. The main 
source of broadening is the increasing band-width for 
higher energy lattices bands (represented by the shad- 
ows around the peaks in Fig. [5]). Due to large lattice 
periodicity in the y direction, the tunneling rate is small 



even for the second excited band corresponding to the 
lowest energy resonance at WLar = Wy = 2ir ■ 55 kHz but 
it increases rapidly for excited bands in the other direc- 
tions (as represented in Fig.EJ ■ Magnetic field and lattice 
depth fluctuations are broadening the resonances, at the 
10 kHz level for data shown in Fig. [5J Another cause 
of broadening is cascading dipolar processes correspond- 
ing to excitations to higher bands, whose energies are 
not equally spaced due to anharmonicity (at the 2 kHz 
level). Cascading is indeed observed in our experiment 
in which states down to \ms = —2) are produced at long 
time (Fig. [T]b). Inhomogeneous broadening of the band 
structure due the finite size of the lattice beams is of the 
order of 500 Hz. 

The lowest frequency resonance has a small broaden- 
ing due to tunneling (100 Hz), allowing to reveal the 
role of on-site interactions. This resonance occurs at 
^Lar = <-<Jy corresponding to an orbital excitation of 2 
units (2hujy) and a change in Zeeman energy of 2gsHBB 
(the coupling operators are "V% and V3). In the rest of 
this letter we focus on this resonance. For these mea- 
surements, shown in Fig. [3][a) after 11.5 ms of relax- 
ation, the magnetic field and the laser intensity were ac- 
tively stabilized to reduce technical broadening. These 
data were also registered at slightly lower lattice depth 
than for Fig. [2] data (corresponding to uj y = 45(1) kHz). 
To explain the observed data, the contact interactions 
U (ri — Y2) = Uq5(yi — T2) have to be taken into account 
( u o = J2s t =o,2A,6 9s M Ps M , p s M = \Sm){S m \ is the 
total spin projector, gs M = ^Trh 2 as M /m and as M the 
scattering length for the different molecular channels). 
The interactions shift the position of the resonances: for 
the lowest one we obtain cjLar = u y H — ^ — - (as sketched 
in Fig. (3Jb)). We calculate the interaction energy shifts 
of the initial and final states and A f using first order 
perturbation theory assuming a harmonic trap. 

For two particles in the same site, the on-site en- 
ergy shift of the initial state \i) is given by Aj = An. 

3 r 

where Ao = (m/h) 2 (ojxOjyOJz)? gg. The final state are 
= |0,2,0)®|6,4) and \f 2 ) = |0, 2, 0}® |4, 4). The spin 
dependent interactions lift the degeneracy and the two fi- 
nal states shifts are A/ 1 = Ao/2and A/2 = (g4/ge) Ao/2. 

As the initial state is more confined than the final ones, 
on-site interactions induce a global red shift of the res- 
onance (about 2.5 kHz) and a splitting of the resonance 
into two peaks (separated by 0.85 kHz, see Fig. [3] (a)) 
corresponding to two molecular potentials. We do not 
resolve the two resonances for two atoms per site, likely 
due to broadening associated with residual magnetic field 
fluctuations and lattice inhomogencities. 

In another set of measurements, we have loaded the 
lattice faster (in less than 5 ms). The line-shape of the 
resonance is then completely modified, as shown Fig. 02 
and presents a multi-resonant structure that we inter- 
pret as due to dipolar relaxation within sites containing 
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three atoms. For atoms in stretched states, the on-site 
interaction shift scales with the number of pairs. How- 
ever, it is not the case when the spin degree of freedom 
is released and the previous result with 2 atoms per site 
cannot be simply scaled by the number of pairs. For this 
specific case with three atoms per site we will now show 
that DDIs couples the initial state to four different fully 
symmetric states, two of them having entangled spin and 
orbit. The observed multi-resonant structure originates 
from these states. 

To solve the case of 3 atoms per site, we have used the 
Jacobi coordinates: the orbital part of the three-body 
wave function can be described by center of mass co- 
ordinate and two relative coordinates (see supplemental 
material) . As in the previous two-body problem, the cen- 
ter of mass motion is decoupled, the total magnetization 
decreases by two units as the orbital energy increases by 
2hu y . The initial state has maximal magnetization M = 
9 and each particle is in the ground state of the harmonic 
oscillator. The shift of this initial stretched state is Aj = 
3A . We find four final bosonic states without center of 
mass excitation, a total magnetization M — 7 and exci- 
tation energy of 2Tiu y . The shifts are given by the expres- 
sions A/i = (9/8) A , A /2 = (9/4) A , A /8 = (o - b)A 
and A/4 = (a + 6)A with a = (150^4 + 147p 6 )/(176g 6 ) 
and b = (5092.9! + 508.g 4 .96 +4201^) 1 / 2 )/(1763 6 ). It can 
be noticed that, when 174 = ge, only two resonances are 
expected. 

Two of these four states are of the form \0) ® \S), 
where \G) and \S) are symmetrical states of the orbital 
and spin parts respectively. The two others are entangled 
states of orbital and spin degrees of freedom. Similar 
states produced by spin changing collisions have also been 
discussed in the context of Fermi spinor gases [3lJ . 

We find a good agreement between our theory and the 
experiment (see Fig. [3]). This shows that we can pro- 
duce and resolve 3-body states with entangled spin and 
orbital degrees of freedom in the lattice. This may also 
enable interesting experiments where correlated tunnel- 
ing of these exotic states could be studied. 

This data is also relevant from the many-body point 
of view. Indeed, Fig. [3] (a) also shows that provided the 
lattice is loaded slowly, no sites with three atoms are pro- 
duced. This is the signature that a squeezed atom distri- 
bution peculiar to the Mott state has been reached. On 
the other hand, for a short loading time sites with three 
atoms are produced. In Fig. [3] (c) we show the resulting 
increase of the resonance width for loading times shorter 
than 5 ms. This adiabatic timescale is comparable to the 
trap period, which suggests that the breakdown of adia- 
baticity for fast loading arises from a time scale that is so 
short that the cloud does not have time to swell to reach 
the ground state of the system. 

In conclusion, we have observed anisotropic dipolar 
resonances that couple the spin and the orbital degrees 
of freedom of bosons in a 3D lattice. We have focused on 



the narrowest resonance which can be used to engineer 
few-body entangled states and to probe the atom number 
distribution in the Mott state. Dipolar resonances pro- 
duce correlated states whose tunneling properties would 
be interesting to study. An interesting prospect would 
also be to reach the coherent regime. For example, start- 
ing from atoms in |— 2) one might expect Rabi oscilla- 
tions coupling states of different magnetization and dif- 
ferent orbital states. Finally, we also stress that away 
from dipolar resonances, the magnetization is fixed, and 
the gas is stable. The system is then a good candidate 
to study lattice magnetism at constant magnetization, 
with interesting connections to quantum magnetism [32| . 
The magnetic stability away from dipolar relaxation res- 
onances displays similarities with the stability of quan- 
tum dots discussed in the context of quantum informa- 
tion computing [33| . 

Acknowledgements: LPL is Unite Mixte (UMR 7538) 
of CNRS and of Universite Paris 13. We acknowledge 
financial support from Conseil Regional d'lle-dc-France, 
Ministcrc dc l'Enseignement Supcrieur et de la Recherche 
(CPER) and IFRAF. A. C. thanks CNRS for financial 
support. We thank J.V. Porto and M. Gajda for their 
suggestions and their comments during the preparation 
of the manuscript. 



[1] M. Baranov, L. Dobrek, K. Goral, L. Santos, and 
M. Lewenstein, Physica Scripta 2002, 74 (2002). 

[2] M. Baranov, Physics Reports 464, 71 (2008). 

[3] T. Lahaye, C. Menotti, L. Santos, M. Lewenstein, and 
T. Pfau, Reports on Progress in Physics 72, 126401 
(2009). 

[4] J. Stuhler, A. Griesmaier, T. Koch, M. Fattori, T. Pfau, 

S. Giovanazzi, P. Pedri, and L. Santos, Phys. Rev. Lett. 

95, 150406 (2005). 
[5] T. Lahaye, J. Metz, B. Frohlich, T. Koch, M. Meister, 

A. Griesmaier, T. Pfau, H. Saito, Y. Kawaguchi, and 

M. Ueda, Phys. Rev. Lett. 101, 080401 (2008). 
[6] G. Bismut, B. Pasquiou, E. Marechal, P. Pedri, 

L. Vernac, O. Gorceix, and B. Laburthe-Tolra, Phys. 

Rev. Lett. 105, 040404 (2010). 
[7] G. Bismut, B. Laburthe-Tolra, E. Marechal, P. Pedri, 

O. Gorceix, and L. Vernac, Phys. Rev. Lett. 109, 155302 

(2012). 

[8] K. Aikawa, A. Frisch, M. Mark, S. Baier, A. Rietzler, 
R. Grimm, and F. Ferlaino, Phys. Rev. Lett. 108, 210401 
(2012). 

[9] M. Lu, N. Q. Burdick, S. H. Youn, and B. L. Lev, Phys. 

Rev. Lett. 107, 190401 (2011). 
[10] M. Lu, N. Q. Burdick, and B. L. Lev, Phys. Rev. Lett. 

108, 215301 (2012). 
[11] A. Auerbach, Springer- Verlag, N.Y p. 255 (1994). 
[12] A. Micheli, G. K. Brennen, and P. Zoller, Nature Phys. 

2, 341 (2006). 

[13] A. V. Gorshkov, S. R. Manmana, G. Chen, J. Ye, 



5 



E. Demler, M. D. Lukin, and A. M. Rey, Phys. Rev. 

Lett. 107, 115301 (2011). 
[14] D. Peter, S. Miiller, S. Wessel, and H. P. Buchler, Phys. 

Rev. Lett. 109, 025303 (2012). 
[15] C. Menotti, C. Trefzger, and M. Lewenstein, Phys. Rev. 

Lett. 98, 235301 (2007). 
[16] K. Goral, L. Santos, and M. Lewenstein, Phys. Rev. Lett. 

88, 170406 (2002). 
[17] E. G. Dalla Torre, E. Berg, and E. Altaian, Phys. Rev. 

Lett. 97, 260401 (2006). 
[18] Y. Kawaguchi, H. Saito, and M. Ueda, Phys. Rev. Lett. 

96, 080405 (2006). 
[19] L. Santos and T. Pfau, Phys. Rev. Lett. 96, 190404 

(2006). 

[20] K. Gawryluk, M. Brewczyk, K. Bongs, and M. Gajda, 

Phys. Rev. Lett. 99, 130401 (2007). 
[21] B. Sun and L. You, Phys. Rev. Lett. 99, 150402 (2007). 
[22] T. Swislocki, T. Sowihski, J. Pietraszewicz, M. Brewczyk, 

M. Lewenstein, J. Zakrzewski, and M. Gajda, Phys. Rev. 

A 83, 063617 (2011). 
[23] B. Pasquiou, G. Bismut, E. Marechal, P. Pedri, 

L. Vernac, O. Gorceix, and B. Laburthe-Tolra, Phys. 

Rev. Lett. 106, 015301 (2011). 
[24] B. Pasquiou, G. Bismut, Q. Beaufils, A. Crubellier, 

E. Marechal, P. Pedri, L. Vernac, O. Gorceix, and 

B. Laburthe-Tolra, Phys. Rev. A 81, 042716 (2010). 
[25] B. Pasquiou, E. Marechal, G. Bismut, P. Pedri, 

L. Vernac, O. Gorceix, and B. Laburthe-Tolra, Phys. 

Rev. Lett. 106, 255303 (2011). 
[26] B. Pasquiou, E. Marechal, L. Vernac, O. Gorceix, and 

B. Laburthe-Tolra, Phys. Rev. Lett. 108, 045307 (2012). 
[27] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and 

P. Zoller, Phys. Rev. Lett. 81, 3108 (1998). 
[28] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and 

I. Bloch, Nature 415, 39 (2002). 
[29] Q. Beaufils, R. Chicireanu, T. Zanon, B. Laburthe-Tolra, 

E. Marechal, L. Vernac, J.-C. Keller, and O. Gorceix, 

Phys. Rev. A 77, 061601 (pages 4) (2008). 
[30] P. L. Gould, G. A. Ruff, and D. E. Pritchard, Phys. Rev. 

Lett. 56, 827 (1986). 
[31] N. Bornemann, P. Hyllus, and L. Santos, Phys. Rev. Lett. 

100, 205302 (2008). 
[32] R. Barnett, D. Petrov, M. Lukin, and E. Demler, Phys. 

Rev. Lett. 96, 190401 (2006). 
[33] M. Paillard, X. Marie, P. Renucci, T. Amand, A. Jbeli, 

and J. M. Gerard, Phys. Rev. Lett. 86, 1634 (2001). 



Supplemental material 

Transition matrix elements. Here we present the re- 
sult of the calculation of the matrix elements related to 
the dipolar relaxation resonances involving two atoms per 
site. The site frequencies arc u> x /2tt = 170 kHz, lo v /2tt = 
55 kHz and ui z /2ir — 100 kHz. Channel 1 corresponds to 
the final spin state |6,5) raol = ^(|3,2) ato + |2,3) ato ). 
Channel 2 and 3 that correspond to final spin states 
6,4) mol and |4, 4) mol respectively are almost degener- 
ated. For data presented in Fig. 2 , the splitting due 
to contact interaction is not resolved. Therefore we 
present in the Table the coupling to the atomic final 

state |2,2) ato = ^\6A) mol - ^ |4,4) mol instead of 
the coupling to the molecular spin states. The final state 
\n x ,n y ,n y ) corresponds to the excitation of the relative 
motion along the three different axis. We give the fre- 
quencies of the resonnance and the couplings (/| Vdd \i) 
where \i) = |0, 0, 0) €5 |6, 6) is the initial state. 



channel 1 



Final 


Resonance 


Coupling 


Coupling 


state 


(kHz) 


B||Ox (Hz) 


B||Oy (Hz) 


10,2,0} 


110 








|0,1,1) 


155 





412 


|0,0,2) 


200 








channel 2 & 3 


Final 


Resonance 


Coupling 


Coupling 


state 


(kHz) 


B||Ox (Hz) 


B||Oy (Hz) 


1 0,2,0) 


55 


107 


47 


|0,1, 1) 


77.5 


238 





|0, 0, 2) 


100 


230 


152 


10,4,0} 


110 


164 


35 


|1,1,0) 


112 








|0,3, 1) 


132.5 


227 





|1,0,1> 


135 





305 


|0,2, 2) 


155 


62 


84 


|0,6, 0) 


165 


190 


28 


1 1,3,0) 


167.5 








|2,0, 0) 


170 


38 


278 


|0,1,3) 


177.5 


208 





|0,5, 1) 


187.5 


223 





1 1,2,1) 


190 





176 



Calculation details for coupled three body states. We 
present here the calculation of the dipolar relaxation res- 
onance spectrum (see Fig. 3) when 3 atoms are loaded in 
the same site. Initially, 3 atoms occupy the ground state 
of an anisotropic harmonic site with frequencies uj x , u> y 
and L) z . DDIs induce magnetization changing collisions 
with excitations to higher orbital states. We consider 
in the problem only orbital excitations along the Oy di- 
rection. Indeed, under our experimental conditions, cu y 
is the smallest frequency and the registered lowest fre- 
quency resonance corresponds to excitation of the 2huj y 
(without central of mass excitation) that comes with a 
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decrease of the total magnetization by two units. We 
describe the three body states by means of the states 
| rii, ri2, ti 3 ; mi, m 2 , m 3 ) where is the Oy harmonic os- 
cillator quantum number and m, the spin projection for 
particles i = 1,2,3. 

The orbital energy is E = %(ni + ?i 2 + n 3 )uj y + Eq and 
total magnetization M = 7711+7712 + 7713. The initial state 
is \i) = |0, 0, 0; 3, 3, 3} with energy Eq. 

The orbital final states with AE = 2huj y are 1 2 , , 0) , 
|0, 2, 0), |0,0,2), |1,1,0), 1,0, 1) and |0, 1,1). With these 
states, we can construct 3-body states without center of 
mass excitations. For this purpose, we use Jacobi co- 



ordinates R 



yi+y-i+y?, 
3 



iv\ - 2/2), P 



_ V1+V2 



?/3 and we construct three operators cir = ai +^| +a3 , 



V2 



and 



( 0.1+0,2 
\ 2 



a-3 



where d, 



~~zh~ ( & rrtrPvi) • nrs t operator acts on the cen- 
ter of mass coordinates, and the two others solely on the 
relative motion. Using these operators, the hamiltonian 

tftf A-tftf\A-En. From the 



now reads H = fou y (a R a R + dj,dj, + a [ p a* p/ 



(a*) 2 (a+) 2 



orbital ground state |0, 0, 0} we can apply —^f and 

alap to generate three orthogonal states without center 
of mass excitations (these three states are a linear com- 
bination of the six final orbital states listed just before). 

The initial spin state is |3, 3,3) with magnetization 
M = 9. The possible final spin states having magne- 
tization 7 are |3,3,1), |3,1,3), |1,3,3), |3,3,2), |3,2,3) 
and |2,3,3). Combining the three orthogonal orbital 



states, we thus obtain 18 orthogonal 3-body states with 
E = 2hw y + Eq, no center of mass excitations and total 
magnetization M = 7. These states describe distinguish- 
able particles. We symmetrize these states in order to 
satisfy bosonic statistics and we finally obtain 4 orthog- 
onal bosonic states without center of mass excitation. 

Effect of contact interactions. These states are de- 
generated without interaction. Spin dependent contact 
interaction fully remove the degeneracy. For three par- 
ticles, it reads U(ri, r 2 ,r 3 ) = Ui 2 S(ri - r 2 ) + ^23^2 - 
r 3 ) + *7i3<K r i - r 3) where E/y = V'; 0,2,4,6 0s p s » 
Psij = (Sij\ is the total spin projector for the par- 
ticles i and j. gs u = 4?T?i 2 as M /m and as M the scat- 
tering length for the different molecular channels. For 
chromium ag = 102.5 ao and 04 = 64 ao [24] where ao is 
the Bohr radius. Within first order perturbation theory, 
we diagonalize analytically the 4x4 Hamiltonian, obtain- 
ing the energy shifts A/ n reported in the paper and the 
final orthogonal states |/„) . Among the four states two of 
them are product state of orbit and spin, the two others 
are entangled states of orbit and spin. With the exper- 
imental parameters corresponding to Fig. 3, we obtain 
(A fl , A f2 , A h ,A U ) = (20.7, 16.8, 10.4, 8.4) kHz 

DDIs couple the initial state to the four final states. 
The matrix elements of the dipolc-dipole interaction 
arc (f n \3Vdd\i) = a n V. We obtain V = 185 Hz and 
{ai,a 2 ,a 3 , a 4 ) = (0.51, 0.43, 0.44, 0.59). 



